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ABSTRACT 


Methods to improve the performance of reflux heat 
pipes for heat recovery applications have been 
examined both analytically and experimentally. 

A survey of literature models for the estimation 
of reflux heat pipe transport capacity was per- 
formed and these various models were compared 
with experimental data. A high transport capacity 
reflux heat pipe was developed that provides up to 
a factor of 10 capacity improvement over conven- 
tional open- tube designs; analytical models were 
developed for this device and incorporated into a 
computer program "HPIPE". Good agreement of the 
model predictions with data for R-ll and benzene 
reflux heat pipes was obtained. 
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1.0 INTRODUCTION 

This report discusses analyses and experiments directed towards the develop- 
ment of an improved reflux heat pipe for heat recovery applications. The 
result of this work has been a "cocurrent flow" heat pipe that has a transport 
capacity of about 10 times that of a conventional open-core heat pipe or two- 
phase thermosyphon. This performance increase is the result of introducing 
a relatively simple plate/tube separator assembly into the heat pipe vapor 
core. This separator system provides cocurrent vapor-liquid flow in all heat 
pipe sections and eliminates the vapor/liquid countercurrent shear that com- 
monly limits reflux heat pipe capacity and causes vibration and thermal in- 
stability. 

Section 2 of this report discusses and compares various literature models of 
reflux heat pipe transport capacity and describes the analytical models 
developed for cocurrent flow heat pipes. The thermal requirements of heat 
pipes for heat recovery applications are outlined and compared to the various 
heat pipe capacities predicted for conventional and cocurrent flew heat pipes. 

Section 3 presents a discussion of experimental tests performed in support of 
the analytical studies and a comparison of these tests with model predictions. 
Section 4 is a summary of the work performed and a discussion of further areas 
meriting .study. 
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2.0 ANALYSIS 

This section describes past and present art in the field of reflux heat pipes 
and discusses the transport requirements of reflux heat pipes used for heat 
recovery applications. Expressions proposed by various investigators for 
estimation of reflux heat pipe transport capacity are presented and compared 
with experimental data. A new reflux heat pipe design is discussed that has 
a particularly high transport capacity and an analytical model for the device 
is presented. 

2.1 Reflux Heat Pipes 

2.1.1 Prior Art 

The application of reflux heat pipes or two-phase thermosyphons to energy- 
related problems in a sense began prior to the development of heat pipes in 
that the concept of using closed heat transfer tubes partially filled with a 
vaporizable fluid was the subject of U. S. patent activity as early as the 
1920' s. Figure 2.1.1 shows the "Heat Transfer Means" presented by F. W. Gay 
in 1929 for recuperator heat recovery. ^ The air inlet and stock exhaust 
ducts were coupled by 

"sealed tubes containing a volatile liquid, said tubes having 
their lower ends in contact with the outgoing hot stack gases 
and their upper ends in contact with ingoing air, so that 
said volatile liquid boils in the lower ends of the tubes 
under the influence of the heat of the stack gases to produce 
a hot vapor which rises to the upper ends of the tubes, the 
heat thus carried to the upper tube ends being radiated there- 
from to the air, and such transfer of heat operating to cool 
and condense said vapor for return to the lower tube ends, so 
that the cycle of heat transfer operations thus produced is 
con t inuous ly repeated." 
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FIGURE 2.1.1 


REFLUX HEAT PIPE OR TWO-PHASE THERMOSYPHON PROPOSED FOR HEAT 
RECOVERY APPLICATIONS IN 1929. 
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( 2 ) 

In a second patent issued in 1930, Gay applied these reflux tubes to the 
cooling of underground electric transmission lines. 

In 1965, E. L. Long patented a refluxing two-phase heat transfer tube ^ 3 ^ 

used for permafrost stabilization which was essentially identical to the heat 

transfer tube of Gay in 1929, while in 1973. E. D. Waters patented the 
(4) 

Cryoanchor a second permafrost stabilization device relying on a reflux- 
ing two-phase liquid/vapor cycle. Other permafrost- type reflux devices are 
referenced by O' Byrne. ^ 

The first mention of capillary wick pumping as a means of liquid distribution 
in a two-phase heat transfer tube was in a 1938 patent ^ and the 1944 patent 
of Gaugler, and during the 1960's and early 1970' s, it was customary to 
identify two-phase heat transfer tubes without a capillary wick and with sub- 
stantial fluid charges as "two-phase thermosyphons" while those relying on 

capillary action and thin-film evaporation and condensation were called "heat 
(7) 

pipes". This distinction has gradually been lost, and numerous references 
can now be found to devices which would have once been called thermosyphons 
that are now called heat pipes. Many of these reflux tubes are substantially 
overloaded but employ wicking or grooves for capillary liquid transfer, and 
hence the distinction is also being lost from an operational standpoint. In 
this report, all these heat transfer devices will be termed "heat pipes". 

R. A. Ares obtained a patent in 1972 on an air-to-air heat pipe heat recovery 

system similar to Gay's that was for the precooling of air for large refrigerated 

(8) 

vaults and produce ripening rooms. The duct-to-duct heat pipes were used 
in combination with a vapor compression cycle refrigeration unit to ostensibly 
reduce the refrigeration load required for make-up air. In 1973, J. E. Runyan 
and G. M. Grover of Q-Dot Corporation patented a circumferentially grooved 
heat pipe for near-horizontal disposition in air-to-air heat recovery heat 

(9) 

exchangers which used a cross- shaped extrusion to separate the vapor and 
liquid and thereby reduce wave action, slugging, and other transport-reducing 
liquid-vapor interactions (Figure 2.1.2). This was followed in 1974 by a 
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second air-to-air neat pipe heat exchanger patent by H. G. Barkmann of Q-Dot 
Corporation,* 10 * in which the heat pipes were again very nearly horizontal, 
and could be regulated in transport capacity by tipping. This patent also 
covered heat pipes with no internal wick or grooving, and was a continuation 
of a patent originally abandoned in 1969. Both of these patents indicated a 
nominal liquid load of about one-third the total heat pipe volume. 

In 1975, R. L. Pessolano et al received a patent on an air-to-air or fluid- 
to-air heat pipe heat exchanger for domestic or industrial heating systems 
that is very similar to Gay’s except that the heat pipes are specified to con- 
tain a capillary material for liquid transport and to not be " dependent 

on gravity for operation.”* 11 * 

M. A. Ruch and G. M. Grover of Q-Dot Corporation presented a paper at the 

1976 International Heat Pipe Conference describing heat pipe thermal recovery 

(12) 

applications and design features of heat pipe heat exchangers. The heat 

pipes described as manufactured by their company were not like those described 
in the patents (9) and (10), but rather were characterized as filled with 
liquid to 50 to 70% of free volume and as having a small-diameter tube 60 to 
80% of the overall length within the larger diameter tube to transport liquid 
back to the evaporator section. Transport capacities with halogenated hydro- 
carbon fluids at 2 to 10 atmospheres were stated to be in the range of 1500 to 
2500 watts for a 2.3 cm I.D. tube, and three to four times greater with the 

return artery. A patent was awarded in 1977 to Grover as a result of this 

work. 

In summary, heat pipes or heat pipe-like two-phase heat transfer tubes have 
been used or proposed in the United States for heat recovery applications 
since 1929 or earlier, but have only received substantial commercial attention 
over the past ten years. Improvements to the basic heat pipe for reflux opera- 
tion have been relatively few, with the use of capillary wick or grooves as a 
liquid transport means being one of the most significant advancements over the 
past 40 years. Some preliminary work has been done on improving transport 
capacity, but as will be discussed in the following section, most analytical 
and experimental work is still centered around modeling of the basic reflux 
tube as described by Gay and others. 
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2.1.2 Transport Models 

The prediction of reflux heat pipe transport capacity has received significant 
attention over the past several years, but there is as yet no satisfactory 
analytical model for this heat transfer device. It has been possible through 
dimensional analysis to establish empirical operating limits, and a summary of 
these models follows. This section closes with a comparison of each model with 
measured transport capacities for the 2.54 cm I.D. open core heat pipe described 
in Section 3. 

A substantial number of the capacity models for vertical reflux operation liken 
the failure of a reflux heat pipe to the phenomena called "flooding" in the 
chemical process industry. In flooding, the falling annular film develops 
surface waves which grow to the point that the entire wave is blown back 
towards the condenser leading to evaporator dry-out. 


Wallis Model 

Wallis et a i d^,14) assujne( j that the critical mass flow was related to a 
balance between gas inertial forces and hydrostatic forces, and that the 
critical axial heat flux density Q for turbulent flow in both phases was 

£lW 

given by 



2 

(W/cm ) 


( 2 . 1 - 1 ) 


where D is the characteristic dimension of the heat pipe evaporator/ad iabatic 
cross section, h f ^ is the latent heat of vaporization, and and are the 
liquid phase and vapor phase densities, respectively. The constant C w is a 
constant that is less than or equal to 1.0, the exact value depending on end 
effects associated with the transition between evaporator and condenser. For 
a large condenser-to-evaporator cross-section ratio and a sharp-edged 
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transition, C - 0.725. For heat pipes in which the wetting angle is close 

to zero and the evaporator and condenser have equal cross-sectional areas, the 
constant may be very close to 1.0. 

Kutateladze Model 

The Wallis model does not account for surface tension effects or viscosity 
effects. A second approach to modeling flooding has used the Kutateladze 
number, which expresses the balance between surface tension, hydrostatic, and 
inertial forces. 


Cy. = K v + (2.1-2) 

where 

K v = p v V Cga(p £- p v ,] ' 1 ; (2.1-3) 

K £ = V[ g °(V P v )] * ; (2.1-4) 

and where j v and are the volumetric flows of vapor and liquid per unit 
cross-sectional area, and a is the surface tension. 


Solving (2.1-2) for the case in which the mass flows of vapor and liquid are 
equal, the critical heat flux density is found to be 


*ak 


i 

h fg c k (ofotyO) 
(■•©")■ 


(2.1-5) 


The theoretical value of is yf 3.2. The Wallis and Kutateladze models of 
critical flux are not independent, since the ratio of (2.1-1) to (2.1-5) is 
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( 2 . 1 - 6 ) 


where the quantity in the square root is the Bond number, i.e., the ratio of 
hydros tatic-to-surf ace- tens ion forces. 


Bo = [D^P^-P^/a] 15 


(2.1-7) 


Composite Models 

(14) (16) 

Wallis and Makkenchery and Tien and Chung recognized the fact that 
Q is a phenomenologically more appropriate limit for large heat pipes where 
enclosure effects disappear while is more appropriate for small tubes 
where the volume per unit contact area is much less. Tien and Chung modified 
C^ for use in equation (2.1-5) as follows to bridge the gap between the two 
models: 


C, t = C k *tanh (0.5 Bo"*) (c^ = J 3.2) (2.1-8) 

In the limit of small Bond numbers, when the modified Kutateladze number C 
is substituted into (2.1-5), the resulting flux limit is of the Wallis form 
[equation (2.1-1)], while for large Bond numbers, C tends towards C^ and 
the critical flux tends towards the Kutateladze limit. 

The modified Kutateladze number C kfc does not represent experimental data in 
the transition region as well as might be desired because of the sluggish 
response of equation (2.1-8) to Bo. A more satisfactory transitional 
Kutateladze number which also provides Q for small Bo when substituted 
into (2.1-5) is 

C ks = C k* t 1-Exp (* 

2 

where a = (c /C, ) . 

w k 


(2.1-9) 
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A comparison of these models with a large body of experimental reflux tube 
and open thermosyphon data is given in Figure 2.1.3. The peak flux data have 
been ratioed with respect to equation (2.1-5), using from equation (2.1-8) 
and equal to 0.725 and 0.9 in equation (2.1-9). The data are represented 
more accurately with the proposed empirical model 


where 


ak 


c 2 h. [p 2 ga(o.-p )] 
ks fg v’ l v' J 


1 ' 


( p / p ,)'] 2 


(w/cm 2 ) 


4 


3.2 [l - EXP(- a 2 Bo)] 


( 2 . 1 - 10 ) 


( 2 . 1 - 11 ) 


However, equation (2.1-10) and the Wallis critical flux model may be of ques- 
tionable value for heat pipes having very small cross sections. It does not 
seem physically correct that the Bond number should become of negligible 
importance for small tubes, since this is the direction of increasing surface 
tension forces. Also, the models do not include any effects of liquid or vapor 
viscosity. A major justification for the use of equation (2.1-10) is that a 
large body of data is correlated quite well by this model. 


Bezrodny et al 

Bezrodny and Beloivan have reported correlations based on extensive reflux 
heat pipe data using water, ethyl alcohol, ad Freon-12. They found two 
flow modes for heat pipes in the 1 to 5 cm I.D. range — below a critical fluid 
fill of about 25% of the evaporator volume, the heat pipe would exhibit one 
heat flux limit, while a slightly larger fill would paradoxically result in a 
20 to 25% lower capacity. This was attributed to a change from annular to 
slug flow at higher mass loadings creating a net decrease in returning conden- 
sate. 

The correlations derived from these data by Bezrodny and Beloivan are based 
on a modified Kutateladze number, where the critical heat flux density Q is 
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TIEI^CHUNG CORRELATION {«qn 2.1.8) 

SUV* EQUATION (2.1.9) 

A DOWTHERM-A, SAKHUJA 
O 4-METHYL PYRDINE, SAKHUJA 
V WATER, KUSUOA AND IMURA 
O ETHYL* ALCOHOL, KUSUOA AND IMURA 
A CARBON TETRACHLORIDE, KUSUOA ANO IMURA 
♦ HEXANE, KUSUOA AND IMURA 
O WATER, COHEN AND BAYLEY 



FIGURE 2.1.3. A COMPARISON OF FLOODING DRY-OUT DATA WITH TWO 
HYDRODYNAMIC MODELS. 
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and 


where 



Cp^(P t -P v )]' 


Pg/CgofPj-Py)]* 5 


( 2 . 1 - 12 ) 

(2.1-13) 

(2.1-14) 


The modified Kutateladze number is dependent on the system saturation pres- 
sure P through a dimensionless pressure k , and falls in the following ranges 

^ 3 *" 4 

depending on the value of K . For 5(l0 ) s k ^ 5(l0 ), 

P P 


5.65/K°* 13 s cj^rs 6.02/K 0 ' 13 


(2.1-15) 


For the higher pressure range 5(l0 4 ) s K^ s 3.5(l0^), it was found that 


1.35SC 


2 

kb 


1.6 


(2.1-16) 


Bezrodny and Sakhatsh performed additional tests on tilted water heat pipes of 
(18) 

1.2 to 3.8 cm I.D. Typical data are shown in Figure 2.1.4. The bistable 

operation discussed earlier at higher tilts is quite obvious. For the tilt 
angle 9 in the range 5 to 30°, the data were correlated with the formula 

Q = 3300 P* e °* 54 S .°' 44 (w/cm 2 ) (2.1-17) 

ac s e 

with vapor pressure P g in MPa and the length in meters. 

Kusuda and Imura 

Kusuda and Imura have developed a model for vertical reflux transport capacity 

that is based on a semi-empirical physical model rather than dimensional 
(19) 

analysis. A related approach is given in reference 19a. This relationship 

based on data gathered at 1 atmosphere vapor pressure, gives the critical heat 
flux density Q for tubular heat pipes as 

clC 




ft cure 2 . 1.4. 


EFFECT f >F TILT FROM HORIZONTAL ON REFLUX HEAT PIPES OF BEZRODNY AND 
SAKHATSH . 
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ac 


0.4 p h, 
M v fg 




(W/cm^) 


(2.1-18) 


where D = heat pipe diameter, u = viscosity, h^ = latent heat per unit mass, 
L g = evaporator length, C = specific heat, and AT g = condensate subcooling 
relative to the vapor core. 


University of Stuttgart 

Abhat, Nguyenchi, Groll, and others at the Institut fur Kernenergetik have been 

studying smooth-wall and screened-wall reflux heat pipes and have published 

(20-23) 

several papers on this work. An important contribution of this work has 

been detailed study of the relationship between axial and radial heat flux 
limits. For short heat pipes operated with low vapor pressure fluids, the heat 
flux at the heat pipe surface may lead to film boiling or depriming prior to 
any axial heat flux limit being exceeded. Based on influence parameter studies 
of water-charged screened wick heat pipes, the following power-law relations 

have been found to describe dependence of limiting heat flux on system 

• u, (23) 
variables . 

+3.49 dependence on heat pipe O.D. 

+0.727 on the fractional liquid fill of the evaporator 

+0.424 dependence on sin 9 

-0.404 dependence on evaporator length 

-1.0 on total fluid volume 


2.1.3 Comparison of Models 

Figures 2.1.5 and 2.1.6 show the various models described in this section com- 
pared with experimental data gathered at this facility with the 2.54 cm I.D. 
reflux heat pipe shown in Figure 3.1.4. For these tests, the heat pipe core 
was open. The measured critical heat fluxes should be somewhat lower than 
the predicted values since heat pipe burnout was identified as the point at 
which the overall heat pipe conductance began to decrease, since the general 
data trend was an increasing overall conductance unless burnout was imminent. 
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25 50 75 100 


TEMPERATURE (°C) 


FIGURE 2.1.5. CRITICAL AXIAL HEAT FLUX PREDICTIONS FOR A 

2.54 CM I.D. HEAT PIPE CHARGED WITH BENZENE. 



CRITICAL HEAT FLUX (worn) 
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TEMPERATURE <°CT 


FIGURE 2.1.6. CRITICAL AXIAL HEAT FLUX PREDICTIONS FOR A 
2.54 CM I.D. HEAT PIPE CHARGED WITH R-ll. 
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In the case of benzene, the burnout was very rapid, with a sharp reduction 
in conductance , while with R-ll, the conductance would often decrease 
gradually so that burnout was difficult to identify. The data agrees well 
with equation (2.1.9) if is assumed to be in the range 0.7 s s 0.95. 

2.1.4 Transport Requirements 

Up to this time, the application of heat pipes to heat recovery has been 
primarily in the area of air-to-air recuperators. Heat from exhausting hot 
gases (such as from furnaces, ovens, dryers, incinerators, gas turbines, stale 
building air, etc.) is reclaimed and used to heat incoming cold air as 
described in Section 2.1. In many of these air-to-air applications where the 
overall temperature difference between the hot and cold streams is small, the 
performance of the heat exchanger is primarily limited by the thermal resistances 
associated with the finned exteriors. As the temperc .ure difference between the 
two streams increases, however, it is possible to exceed either the internal 
surface vaporization limit of the heat pipe or its axial transport capacity. 

These limits will also be approached with greater frequency when the heat pipe 
is used in situations where exterior heat transfer coefficients are larger, such 
as in air-to-f luid, liquid- to-air, or liquid-to-liquid heat recovery applica- 
tions. 

Various correlations exist for heat transfer coefficients associated with air- 

to-air finned tube heat exchangers. ^ 24 26 ^ An empirical model for triangular- 

pitch heat exchangers that has been found to be reasonably accurate is that of 
(24) 

Briggs and Young, 


h 


0.1378 Re°' 718 Pr 1/3 K (S/L)°' 296 

a 


D 

r 


(W/cm K) 


where Re = D G /u . 

r m a 


(2.1-19) 
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The Reynolds number P.a is based on the O.D. of the tube, D^, the mass flow 

rate per unit area at minimum heat exchanger cross section, G , and the gas 

m 

viscosity, u . The heat transfer coefficient is, in addition, dependent on the 

& 

Prandtl number, Pr, the air gap between fins, S, the fin length, L, and the 
air conductivity, K . The exterior conductance per unit length of a 2.54 cm 
O.D. heat pipe with L = 1.27 cm is shown in Figure 2.1.7 for fin densities 
of 2.36 cm * and 4.72 cm * (neglecting heat transfer to the tube itself). For 
each fin density, a range of conductances is shown that corresponds to a fin 
thickness span of 0 to 0.10 cm. Table 2.1.1 summarizes the radial heat flux 
densities and axial heat fluxes per unit evaporator length for typical flue 
gas velocities and various gas-metal temperature differences. The lower end 
of the temperature difference range is indicative of air conditioning applica- 
tions, the middle temperature differences are indicative of industrial air-to- 
air recuperators, while the upper end is indicative of air-to-liquid systems 
for the preheating of process fluids by flue gases. 

Typical evaporative heat flux densities for low temperature difference applica- 

2 

tions range from less than 1 to about 2 W/cm . Axial heat flux requirements 

range from several hundred to less than 2000 watts per meter of evaporator 

length. Industrial air-to-air recuperator-type applications show evaporative 

2 

heat flux densities in the range of 2 to 15 W/cm and axial heat fluxes of a 

few thousand watts to more than 10 kilowatts per meter of evaporator. Air-^to- 

2 

liquid applications show heat flux densities from about 5 to 25 watts/cm and 
axial fluxes from several thousand watts/m to almost 20 kW/m. 

Most of Dowtherms, organic, and halocarbon working fluids show axial reflux 
transport limits for a 2.54 cm I.D. that are in the range of 1000 to 2500 
watts. It is therefore difficult to apply conventional reflux heat pipes to 
heat recovery applications where the driving temperature differences are large, 
such as stack gas recuperators, and/or where the evaporators are longer than 
about 1 meter. A similar or worse state of affairs is found for liquid/liquid 
or liquid/ two-phase process streams where the heat transfer coefficients are 
substantially larger than those found with forced-air and finned-tube systems. 




MAX. FLUE GAS VELOCITY (CM/ SEC) 


FIGURE 2.1.7. CONDUCTANCE PER UNIT TUBE LENGTH FOR A FINNED 
TUBE IN CROSS FLOW FOR FIN DENSITIES OF 
2.36 CM-1 AND 4.72 CM - *. 
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TABLE 2.1.1. HEAT FLUX REQUIREMENTS FOR AIR- 

TO-AIR HEAT RECOVERY HEAT PIPES (D 


Maximum 

Gas 

Velocity 

(cm/sec) 

Air-to-Fin 

Temperature 

Difference 

(K) 

Radial Flux 

2 

Density (W/cm ) 

Per Unit Length 
Axial Heat Flux (W/M) 

2.36 fins /cm 4.72 fins/cm 

2.36 fins/cm 

4.72 fins/cm 

200 

19 

0.4 

0.6 

300 

500 


111 

2.3 

3.6 

1,800 

2,900 


167 

3.4 

5.4 

2,700 

4,300 


278 

5.7 

9.0 

4,600 

7,100 

500 

19 

0.8 

1.2 

600 

900 


111 

4.4 

6.9 

3,500 

5,500 


167 

6.7 

10.4 

5,300 

8,300 


278 

11.1 

17.3 

8,800 

13,800 

750 

19 

1.0 

1.6 

800 

1,300 


111 

5.9 

9.2 

4,700 

7,400 


167 

8.9 

13.9 

7,100 

11,100 


278 

14.8 

23.1 

11,800 

18,500 

(1) Heat 

pipe O.D. = 2 

.54 cm; fin 

length L = 1.27 cm; 



Average air film temperature T = 232 °C; 


Fin width * 0.05 cm; fin effectiveness = 100%. 
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From an economic standpoint, it is further clear that to achieve a given level 
of overall system conductance for a heat pipe heat exchanger, it is more desir- 
able to increase the length of each heat pipe than to install more heat pipes, 
but this direction of economic optimization is restricted by the axial 
transport limit. To expand the use of heat pipes to more varied heat recovery 
applications, it is necessary to develop methods of maximizing heat pipe radial 
and axial heat flux limits and evaporative and condensation heat transfer coef- 
ficients. 

2. 2 Cocurrent Flow Heat Pipe 

As discussed in Section 2.1, the primary factor limiting transport capacity in 
conventional gravity reflux heat pipes is the countercurrent vapor-liquid shear 
process that results in liquid starvation of the evaporator and substantial 
flooding of the condenser surface. To reach heat flux levels beyond this 
capacity limit, a class of heat pipes has been investigated and modeled that 
operate completely in cocurrent two-phase flow. This is accomplished by forming 
what is essentially an optimized "thermosyphon boiler" within the heat pipe. 

The general configuration studied is shown in Figure 2.2.1. 

In the evaporator and adiabatic, condensate is separated from the turbulent two- 
phase mixture produced by boiling by a sheath of impermeable material. For 
most purposes, this sheath is simply a thin-walled tube. The tube is open at 
the lower end and may have an optional slot or set of bleed holes along the 
length in the evaporator to (a) more uniformly wet the evaporator I.D. and 
(b) reduce the overall pressure drop attributable to condensate return. Liquid 
in the condensate tube does not generally boil except at very low heat flux 
levels, and for all practical purposes, it can be assumed that the condensate 
subcooling is more than adequate to forestall a condensate flow reversal. 

In the condenser, the heat pipe cross section has been divided into halves by 
a tight-fitting separator plate that extends to within a few diameters of 
the blind end of the condenser. At the adiabatic-condenser interface, the 
separator plate and condensate return tube are coupled together so as to direct 
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CONDENSER 


ADIABATIC AND 
EVAPORATOR 


± 


FIGURE 2.2. 




SECTION A-A 


EVAPORATOR 
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and transport condensate back to the evaporator from one of the condenser 
half- sect ions. The condensing two-phase flow executes a hairpin- shaped flow 
pattern in the condenser that s.iears stagnant condensate from the condensing 
surfaces and sweeps it towards the return tube interface. The condenser sec- 
tion that is joined hydraulically to the condensate tube partially fills with 
liquid and serves as a reservoir or hydraulic capacitor, reducing the slugging 
and oscillatory behavior often seen in open-section reflux boilers, while the 
high velocity flow created by the separator enhances the condensing heat trans- 
fer coefficient by reducing c ndensate film thickness. The high velocity co- 
current flow also minimizes liquid puddling and reduces the liquid load neces- 
sary to sustain a given mass transport rate. 

A benefit deriving from this type of reflux heat pipe, in addition to high 
heat flux rates, is a substantial body of relevant analytical and empirical 
studies of cocurrent flow that are of use in modeling and optimization. 

2.2.1 Theory: Cocurrent Flow Heat Pipe 

Ope-ation of a heat pipe configured as in Figure 2.2.1 has been modeled pri- 
marily on the two-phase pressure drop expressions developed by Lockhart and 
(27) (28) 

Martinelli and by Owens. The transport capacity is found by determin- 

ing the mass flow rate at which the available hydrostatic head difference 
between the "up" leg and the "down" leg of the flow circuit is balanced by the 
frictional and acceleration pressure drops associated with the two legs, 

AP = AF + AP £ , (2.2-1) 

g a f 

where AP^ is the available hydrostatic head, AP a is the pressure drop resulting 
from fluid acceleration, and AP^ is the pressure drop resulting from friction. 
The "up" leg of the flow circuit in Figure 2.2.1 includes the evaporator and 
adiabatic core sections surrounding the condensate tube and the half-section 
of the condenser with upward two-phase flow. The "down" leg includes the con- 
denser half-section in which the working fluid flow is downward and the inner 
core of the condensate return tube. 
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The variou pressure drops associated with this cocurrent flow heat pipe are 
diagrammed in Figure 2.2.2 and are as follows: 

Evaporator (upflow) 

• frictional pressure drop, R gf 

• acceleration pressure drop, R 

GcL 

Adiabatic (upflow) 

• frictional pressure drop — prior to separator plate, R af 

• frictional pressure drop — with separator plate, R asf 

• constriction pressure drop — with separator plate, R 

a SC 

Condenser (upflow and downflow) 

• frictional pressure drops, R cuf and R c df 

• deceleration pressure differences, R and R 

cua cda 

• expansion and constriction pressure drops at the condenser 
end of separator plate, R cue and R cdc 

Condensate Return Tube (downflow) 

• frictional pressure drop, R tf 

• expansion and constriction pressure drops, R , R 

The following sections outline the methods used for calculating each flow 
resistance. 


In the subsequent analyses, all or most of the flow cross sections are assumed 
noncircular, and it is necessary' to use the hydraulic radius approximation for 
the calculation of Reynolds numbers and friction pressure drops. The expressions 
below are used extensively in the hydrodynamic modeling of heat pipe function 
that follows. 



h u ~ uP 


R. = A /P 
n p w 
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where R^ is the hydraulic radius, is the wetted perimeter, is the per- 
pendicular flow area, m is the mass flow rate, v is the mean flow velocity, 

flV 

p is the fluid density, u is the fluid viscosity, and F is the Fanning fric- 
tion factor, which is a function of Re^. 

For Reynolds numbers less than 2000, the friction factor is given by 

F - 16/Re, n (2.2-3) 

For Reynolds numbers in the range 2000 < Re < 2.5 (10^), the friction factor expres- 
sion used is for hydraulically smooth tube, and agrees with Prandtl's universal 
(29) 

law of friction to within ±0.5%, 


0 21 96 2 ^ — 3 _ ~i 

F = ' . - - 0.261297 (10 ) In (Re w ) + 5.95156 (10 ) 

„ 0.43 h 

Re. 


(2.2-4) 


2. 2. 1.1 Evaporator 


For the evaporator, the pressure drops R , and R were calculated using the 

ei/ -i a \ ea 


homogeneous two-phase flow model of Owens. 


(28) 


In reference 28, Owens presents an analytical model for the pressure drop in 

homogeneous two-phase flow that compares favorably with experimental data for 

air-water mixtures from 14.7 psia to 3000 psia with qualities of about 10 to 

100%, and for vaporizing flows of water ranging from 64 to 381 psia with 

qualities up to 57 percent. If the quality "x" of the vaporizing flow changes 

linearly with axial evaporator position, and if both the liquid and vapor 

densities and dP v /dP are essentially constant, then the evaporator pressure 

drop R , + R is 
ef ea 
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2 2 

where A * p./p - lj B ■ G d(l/p )/dP; C » G /p fl , where l is the evaporator 
K, v e v e x. e 

length, G is the total mass flow per unit cross-sectional area; x and x, are 
e o 1 

the mixture qualities at the evaporator entrance and exit, respectively, and 
(dP/dz) is the pressure gradient at the blind end of the evaporator, z ■ 0, 
with the flow assumed to be completely liquid, i.e., zero quality. 

The pressure gradient term (dP/dz) is not calculated directly, but is instead 
calculated from the pressure gradient at the evaporator exit. This insures a 
consistent pressure gradient transition at the evaporator-adiabatic interface. 
This gradient term is calculated as follows. 


Owens shows that the friction factor for two-phase flow is very nearly equal to 
the friction factor that would be used if the mixture were all in the liquid 
state, having zero quality. Hence, the frictional pressure drop in two-phase 


flow, - (dP/dz) £, is approximately the product of an arithmetic constant and 
a term involving the mean density. 



where F^ is the liquid-state friction factor, is the average local mixture 
quality, and is the evaporator hydraulic radius. Equation (2.2-6) would be 
exact if were replaced by the unknown two-phase friction factor, F fc . The 
best constant for (2.2- 7 ) to achieve continuities at the evaporator/adiabatic 
interface is found by calculating the frictional pressure gradient at the 
evaporator exit using the Lockhart-Martinelli correlation. 
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2. 2. 1.2 Adiabatic 

Two-phase pressure drops in the adiabatic and separator sections have been 

(27) 

modeled using the empirical correlation of Lockhart and Martinelli. 

With this method of estimating two-phase pressure drop, the pressure drop is 
given by either a vapor phase pressure drop multiplied by a correction factor 

AP fcp = R af = Y 2 • (dP/dz) vs • * a (dynes/cm 2 ) (2.2-9) 

where l is the length of the adiabatic section, AP is the two-phase pressure 

drop, (dP/dz) is the pressure gradient if only the gas were flowing through 

vs 2 

the heat pipe cross section, and Y^ is a correction factor dependent on the 
ratio of the single-phase liquid to single-phase vapor pressure gradient. 


Y 2 = f(x; 
v 


where 



(dP/dz) 4s 


/(dP/dz) 

/ vs 


( 2 . 2 - 10 ) 


2 

In the term X , (dP/dz ) ^ is analogous to (dP/dz) in that it is the pressure 
gradient that would result if only the liquid phase were flowing through the 
heat pipe. 


The value of the parameter Y^ also depends on whether the single-component 

partial flows are laminar or viscous on the basis of the respective Reynolds 

numbers. If both partial flows have Reynolds numbers greater than 2000, the 
2 

function Y v is given to within i5% over the range 0.01 < X< 100 by 



(turbulent vapor, 
turbulent liquid) (2.2-11) 


2 

For liquid-phase Reynolds numbers less than 2000, the function Y v is given over 
the same range to within ±8% by 
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v 


/ \0. 218215 / \ 

li + i.i25x ) • fi + yrj 


3.31326 


(turbulent vapor, 
viscous liquid) ^ 2 _x2) 


The vapor phase is assumed to be in turbulent flow at all times. 


The frictional pressure drop R as ff which must be considered when the separator 
plate extends into the adiabatic, is calculated using a procedure identical to 
that above. 


The constriction pressure drop associated with the transition from the adiaba- 

(30) 

trc to the half-section is calculated assuming an abrupt transition. 


asc 


/g 2 (i + ax)\ ( A\ 

= ^ — J . 0.45 ^1 - (dynes/cm 


(2.2-13) 


where G s is the mass flow rate per unit area of separator cross-sectional area, 
and are the flow cross sections in the adiabatic and separator sections 
respectively, and "A" is given in the nomenclature for equation (2.2-5). 

2.2.1. 3 Condenser 

For all design cases, it is assumed that the separator plate ends at the con- 
denser entrance or penetrates into the condenser for some finite fraction of 
the condenser length, and that the condensation process occurs uniformly 
throughout the condenser. As a result of visual observations of heat pipe 
operation using a glass adiabatic section, the separator plate is found not to 
be completely successful in separating the countercurrent flows within the 
condenser. The upward flowing two-phase mixture can entrain substantial amounts 
of condensate through the two slots formed by a loose-fitting separator plate. 

To account for this leakage, the condenser inlet quality is a user-modifiable 
variable. The exit quality X<. is also an input variable, although under normal 
conditions of fluid charge, X^ can be assumed equal to zero. 
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The frictional and acceleration pressure changes in the condenser arc cal- 
culated using the same homogeneous flow model as in the evaporator. However, 
because of mixture deceleration associated with condensation, there is some 
pressure recovery to compensate for frictional loss. This partial pressure 
recovery will always occur in a zone with decreasing quality. 

With a uniform condensation rate, the quality of the mixture in the condenser 
is assumed to change linearly with axial position with the axial rate-of-change 
doubling once the vapor-liquid mixture exits the first condenser half-section. 
Therefore, the two-phase qualities on opposite sides of the separator plate 
are generally substantially different, even where the upward- and downward- 
moving flows meet at the separator plate exit. These qualities at the separa- 
tor plate exit, and X^, are given by 

2 

h - x 2 + 0t 5 - X 2 > jf- (2.2-14) 

c 

l 

and X 4 - X 5 - (X 5 - X 2 ) (2.2-15) 

c 

where 2 c is the penetration depth of the separator plate in the condenser, 

and 2 is the condenser length, 
c 

The pressure drop terms R cu ^» R CU a' R cdf' an( * R cda' are ca ^ cu ^- ate< ^ using equa- 
tion (2.2-5) with the two-phase qualities, section length, and initial pressure 
drops characteristic of the two condenser half-sections. 

The outlet and inlet constriction resistances R and R . associated with 

cue cdc 

the termination of the separator plate in the condenser have been each approxi- 
mated as the pressure drop associated with a 90° mitered corner. 
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R 


cue 



(dynes/cm 2 ) (2.2-16) 


R . 
cdc 



(dynes/cm 2 ) (2.2-17) 


where G c is the mass flow per unit area on either side of the condenser separa- 
tor plate, and F ^ is the friction factor for flow on either side of the separa- 
tor plate, assuming all mass flow to be in the liquid state, 

F . = F (4 m /u a P ) (2.2-18) 

C? \ C X. wc / 

where m^ is the condenser total mass flow rate, and P^ c is the wetted perimeter 
on either side of the separator plate. 

2. 2. 1.4 Condensate Tube 

The frictional pressure drop in the condensate return tube, R^, is calculated 
using equation (2.2-2), 


R 


tf 



(dynes/ cm 2 ) (2.2-19) 


where F fc is the condensate tube friction factor, R^ is the tube hydraulic radius, 
and l is the condensate tube length. 


The constriction resistance R 

tc 

given by 


at the condensate tube-separator interface is 



( 2 . 2 - 20 ) 
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where G fc is the mass flow per unit area in the condensate tube, A fc is the 

condensate tube flow cross section, and A is the flow cross-sectional area 

x 

on one side of the separator. 

The expansion resistance R fce represents the pressure difference arising from 
flow exiting the condensate tube of area A q into the total evaporator cross 
section and thereafter being diverted 180° and constricted down to a smaller 
area A^ that is equal to the cross-sectional area of the evaporator minus the 
tube cross-sectional area. This multiple effect flow resistance has been 
approximated as the sum of 


( 1 ) 

( 2 ) 

(3) 


An expansion from area A q to area (A^ + A^). 

Two 90° mitre bend of cross-sectional area (A + A, ) . 

o 1 

A constriction from area (A + A, ) to area A, . 

o 1 1 


Employing expansion coefficients in references 30 and 31, the pressure drop is 


te 


G 

Ip 


2 

k ( 


1 + 120 B F 


(1 + B) 


0.45 B 3 \ 

1 + B J 


( 2 . 2 - 21 ) 


where G^_ is the mass flow per unit area in the condensate tube, B = A /'A, , 
t o 1 

and F is the friction factor based on the Reynolds number Re = 4m/[u.(p - P )], 

e x. wl wo 

where m is the total mass flow, u is the liquid viscosity and P , and P are 

wl wo 

the wetted perimeters of the evaporator and condensate tube, respectively. 


2. 2. 1.5 Hydrostatic Pressure Difference 

The driving pressure difference for fluid circulation is the difference in 
hydrostatic head between the "up" leg of the heat pipe and the "down" leg. The 
various gravity pressure differences are shown in Figure 2.2.1 as H^. The 
gravity pressure differences for the evaporator and condenser sections of the 
heat pipe are calculated using the expression derived by Owens for a homogeneous 
mixture. This expression is, for the evaporator gravity head H g , 
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gp S, £ e Sin9 


e " (B - A)(X - x ) 
i o 


[•• (^) * (^)] 


( 2 . 2 - 22 ) 


where B and A have been previously defined for equation (2.2.5), and 0 is che 

angle of the heat pipe axis with respect to horizontal. For the condenser 

heads H and H ., I will be replaced ir equation (2.2-22) by the length of 
cu ca e 

the separator plate in the condenser, l , G will be replaced in "B" by G , 

SC 6 C 

and the entrance and exit qualities X ,X, will be replaced with either X_,X, 

o 1 z 3 

or X . ,x c • 

4 5 

In the adiabatic zone and the vapor-phase side of the separator section in the 

adiabatic, the gravity head H is given by 

3 

H = p gZ sin6/(l + AX ) (2.2-23) 

3L )L Z ±. 

The gravity head for the condensate tube is 


H t = P 2, g£ t Sin6 (2.2-24) 

In these expressions, l and £ are, respectively, the lengths of the adiabatic 

3 t 

and the condensate tube. If the condensate reservoir shown in Figure 2.2.1 fills 
up substantially, 2.^ must be corrected for this additional driving force. 

2. 2. 1.6 Conclusion - Hydrodynamic Model 

The overall driving gravity pressure difference is, therefore, 

2 

AP = H+H.-(h + H + H ) (dynes/cm ) (2.2-25) 

g t cd e a cu 


and the sum of the frictional and acceleration pressure difference terms is 
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AP + AP, 

3l £ 


R _ + R + R, + R_ + R + R , + R + 
ef ea af asf a sc cuf cua 


R + R . + R + R ,. + + R + R 

cue cdc cda cdf tc tf te 


(dynes/cm ) (2.2-26) 


A hydrodynaxnically balanced mass flux results in the equality 


AP = AP + AP, 
g a f 


(2.2-27) 


and this condition signifies a valid heat pipe operating point. 

2.2.2 Parametric Analysis 

The cocurrent flow heat pipe model developed in the previous sections has two 
parameters, X^ and , which must be selected on an empirical basis. The 
parameter X^ is the quality of the two-phase mixture exiting the evaporator 
while X^ is the quality of the t .•o-phase mixture entering the condenser. On 
a theoretical basis, X 2 must equal X^ , but, in fact, it is desirable to main- 
tain X 0 independently adjustable to correct for leakage across the condenser 
separator plate. This plate cannot be expected to hermetically seal the 
opposing half-sections and, in fact, experiments in which a glass heat pipe 
section was used showed that considerable leakage does occur. These experiments 
also showed that the quality of the mixture exiting the evaporator approaches 
90-100% near burnout, so that if burnout points are the only operating points 
of interest, the exit qua! ity X^ does not have to be substantially varied. 


Figure 2.2.3 presents predicted transport capacities at 25°C for the "Phase 7" 
heat pipe described in Section 3 using R-ll as working fluid and X^ =0.9 and 
^2 = +0.67. The capacity is shown as a function of the ratio of the condensate 
return tube hydraulic radius R^ to the hydraulic radius of the evaporator two- 
phase boiling section, R^. The figure shows that the optimum ratio of R^ 
to R^ v is not a function of tilt angle and is approximately R^/R^ v = 0.65. 
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FIGURE 2.2.3. CAl.CUI.ATED TRANSPORT CAPACITY FOR PHASE 7 COCURRENT FLOW HEAT PIPE 
AT 10 AND 90° REFLUX TILTS AS A FUNCTION OF THE LIQUID/VAPOR 
HYDRAULIC RADIUS RATIO FOR THE EVAPORATOR SECTION. 
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This optimum ratio will change with changes in fluid properties and lengths 
of each heat pipe section, but it can be expected that in most cases, the 
ratio will be less than or equal to 1.0. 

Figures 2.2.4 and 2.2.5 show the influence of the qualities X^ and X^ on 
predicted performance. Also shown in these figures are the predicted per- 
formances of bare tube thermosyphons of equal diameter. The predicted criti- 
cal heat fluxes at the best estimate values of evaporator exit and condenser 
inlet qualities (X^ =* 0.9, X 2 - 0.67) are bounded within approximately 20% 

by the conservative estimates of X^ = 0.95, X^ = 0.9, and X^ = 0.8, X^ = 0.5. 

It can, therefore, be concluded that the effects of X^ and X^ are measurable, 
but not dramatically so. It should also be noted that the effects of vapor 
qualities on predicted critical heat flux are much less pronounced at lower 
temperatures (lower saturation pressures) . 

2.2.3 Comparison of Fluids for Heat Pipe Applications 

Table 2.2.1 presents the predicted critical heat fluxes for various fluids 

using either a conventional open-core reflux heat pipe or a cocurrent flow 

reflux pipe. In 11 cases a vertical, 2.51 era I.D. tube was assumed. For the 

cocurrent calculations, a separator design, referred to as Phase 10 in Section 

3 of this report, was assumed. Operating temperatures were selected such 

6 2 

that the saturation pressure of the fluid was 2 atm (2.03 x 10 dynes/cm ). 

By performing all the calculations at the same saturation pressure, the vapor 
energy density (vapor density times heat of vaporization) was approximated equal 
for all the fluids examined. 

The fluids examined covered a range from cryogenic liquids to molten metals 
and included commun refrigerants, heat transfer fluids, hydrocarbons, and 
wa.er. Additional fluids of interest may be found in a report by Saaski and 

- u- < 32) 

Owzarski. 
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TEMPERATURE (°G 

FIGURE 2.2.4. THE EFFECT OF EVAPORATOR EXIT QUALITY, X £V , AND 

CONDENSER ENTRANCE QUALITY, X , ON THE PERFORMANCE 
OF A COCURRENT FLOW HEAT PIPE C CHARGED WITH BENZENE. 



CRITICAL HEAT FLUX TKW) 



TEMPERATURE (°C) 

FIGURE 2.2.5. THE EFFE2T OF EVAPORATOR EXIT QUALITY, X ev , AND 

CONDENSER ENTRANCE QUALITY, X c , ON THE PERFORMANCE 
OF A COCURRENT FLOW HEAT PIPE CHARTED WITH R-ll. 



TABLE 2.2.1. CRITICAL HEAT FLUX FOR VARIOUS PLUIDS AT P SAT * 2 ATM 





p V h fg 

p J. h fg 

Wallis* 


Cocurrent* 


Temperature 


/Joules! 

/Joules \ 

Eqn. 2.1-1 

Eqn. 2.1-10* 

Model 

Liquid 

(°C) 

Bo 

\ gm / 

\ gm / 

(Watts) 

(Watts) 

(Watts) 


N 2 

-190 

26.4 

2.32 

148 

1336 

1,131 

14,124 

Ethane 

-75 

15.4 

1.81 

247 

1678 

1,517 

15,762 

hh 3 

-20 

11.6 

2.10 

882 

3836 

3,548 

32,360 

R-ll 

45 

23.9 

2.20 

297 

2025 

1,736 

20,370 

Benzene 

105 

16.5 

1.95 

297 

1929 

1,734 

17,618 

H.O 

Z 

120 

10.7 

2.42 

1952 

6480 

6,041 

52,349 

TiC 1 „ 
4 

165 

22.8 

1.90 

257 

1753 

1,514 

17,045 

Dowtherm-E 

210 

21.1 

2.09 

278 

1904 

1,663 

18,221 

bowtherm-A 

290 

19.7 

2.02 

2 32 

1680 

1,478 

16,124 

Se 

750 

21.9 

2.10 

3533 

8568 

7,446 

66,606 

K 

840 

8. 2 

1 .96 

1233 

4532 

4,286 

34,884 

-Eqns. 2.1-1 and 

2.1-10 are 

based on 

2.54 cm 

bare tube. 




-Cocurrent model 

is based on 

2.54 cm I 

.D. with 

Phase 10 

separator 




* 

Values expressed to units digit for comparison only. 
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As can be seen from the table, liquid energy density, liquid density times 
heat of vaporization, is a dominant factor in determining the performance of 
both the bare tube and the cocurrent heat pipes. This is to be expected 
since at a constant vapor energy density, the capacity of all these devices 
is limited by the rate at which liquid can be returned to the evaporator. 

Higher liquid densities provide greater driving forces to return the liquid 
and higher heats of vaporization yield larger heat transfer rates per volume 
of liquid evaporated. It is interesting to note that for the cocurrent heat 
pipe with the Phase 10 separator, the critical heat flux correlates well 
with the square root of the liquid energy density, [Q(kw)ss34.4 \f (kj/gm)]. 

Titanium tetrachloride is of particular interest in that it is thermally stable 
and relatively noncorrosive in a temperature range of interest for commercial 
heat recovery applications. Tests at Sigma Research, Inc. have shown over 
20,000 hours of operation of TiCl^ in heat pipes made of carbon steel at 
temperatures in the range of 150-170°C. 

Of the non-liquid metals, water shows the highest heat flux capability for both 
open core heat pipes and cocurrent flow heat pipes, with calculated transport 
capacities respectively of 6500-11,000 watts and 52,000 watts. The poorest 
fluid is liquid nitrogen at 1300-1600 watts and 14,000 watts capacity, respec- 
tively, for conventional and cocurrent flow 2.54 cm diameter heat pipes. Most 
refrigerants and heat transfer fluids indicate transport capacities of 1100-2500 
watts for a conventional reflux heat pipe, and 15,000 to 20,000 watts in a co- 
current heat pipe. In general, it can be concluded that a factor of 5 to 10 
improvement in capacity is obtained by eliminating countercurrent vapor/liquid 
shear . 
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3.0 EXPERIMENTAL TESTS 

3.1 Test Apparatus 

Heat pipes developed for this project were tested in one of two test stands 
constructed for this purpose. These test stands consisted of a heat source, 
cooling water supply, means for adjusting the tilt of the heat pipe, and ap- 
propriately placed thermocouples . The first test stand developed used low- 
pressure steam to heat the evaporator. In later tests, a direct coupled water 
loop was used. 

3.1.1 Low-Pressure Steam System 

In this system, shown in Figure 3-1.1, heat was transferred to the heat pipe 
using a water boiler and steam jacket. This method of heating approximated 
a constant-temperature boundary condition and allowed the heat pipe to operate 
in a partially burned out evaporator condition. A previous test system using 
a resistance-heated evaporator was discarded because of hot spot problems 
resulting from local dry-out, and because a constant flux boundary condition 
does not truly reflect the performance of a heat recovery heat pipe. In a 
heat recovery system, the heat pipe wall does not rise indefinitely in tempera- 
ture if a partial dry-out occurs because (1) the gas stream-to-f in temperature 
differential falls, reducing local heat input; (2) the heat pipe wall and 
exterior fins conduct heat circumferentially to areas that are wetted; and 
(3) the available driving temperature differential may not be particularly large 
in many applications. Furthermore, extensive reflux heat pipe test programs in 
the USSR^ 7 ^ have used resistance heating and the adverse effects of tilt on 

capacity discussed in these reports do not coincide with the observed behavior 

( 12 ) 

of air-to-air heat recovery heat pipes produced by U. S. manufacturers, 
while the data collected in this program do agree with the latter report in 
terms of tilt effects on transport capacity. However, it should be recognized 
that a constant- temperature boundary condition is still only an approximation 
of diverse possible operating conditions, and some differences in performance 
can be expected between this simulation and the actual article near burnout. 
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FIGURE 3.1.1. INITIAL TEST SET-UP FOR REFLUX HEAT PIPE TESTS. 
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Heat input to the heat pipe was measured in two ways. The steam boiler supply- 
ing the heat pipe was heated by a water loop, and heat input to the boiler by 
this liquid loop was measured by calorimetric means (mC AT), while the heat pipe 
condenser was cooled by an annular water jacket that also could be used as a 
flow calorimeter at higher power levels. Both sets of thermocouples were 
selected on the basis of identical readings at ice bath temperature and 100°C, 
when used in conjunction with a 0.1°C resolution Doric DS-350 digital thermo- 
couple readout. The Chromel-Alumel thermocouples were directly soldered to the 
inlet and outlet water connections and were well insulated. 

External heat loss from the insulated steam jacket was measured with the 
heat pipe in a dry condition, and was found to be described by a constant con- 
ductance of 1.41 W/ C c. Hence, at a typical water temperature of 65°C, about 
63 watts were lost to a 20°C ambient. 

Prior to testing, the heat pipe was mounted in a large rigid frame and adjusted 
to a vertical orientation using a 61 cm long spirit level. Other angles were 
calculated by direct measurement of differential height from end-to-end of the 
heat pipe. 

3.1.2 Water Loop 

In order to operate at the higher powers necessary to test the high-capacity 
heat pipes developed during this program, a water loop heating system was con- 
structed. This loop is shown schematically in Figure 3.1.2. Also shown in 
this drawing are the locations of the thermocouples used to measure the per- 
formance of the heat pipe. The specific thermocouple numbers given in Figure 
3.1.2 are those used in the Mod-6 through Mod-10 tests; other numbering schemes 
were used in earlier tests. 

A mixinr ”essel, with a volume of approximately two liters, was used to provide 
a free surface. This ensured that the system could not become pressurized. 

The vessel drained into the inlet of a small pump. The capacity of the pump 
was approximately 350 to 425 cubic centimetars per second, depend. ng on the 
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piping arrangement used for the test. Measurements showed the flow rate, for 
any particular arrangement, to be constant over the temperature ranges 
experienced during operation (25-100°C) . 

The discharge of the pump was connected to a heater manifold. This manifold 
initially contained approximately 3000 watts of electrical heaters. As heat 
pipe capacities increased, additional elements were added to provide 6000, 

12,000 and finally 15,000 watts of heat input. The heating elements used were 
commercially available units for hot water heating. 

The mixing vessel, pump, and heater manifold were mounted on a chassis along 
with the necessary electrical controls. This assembly was connected to the 
evaporator heating jacket with foam- insulated high-pressure high- temperature 
hose. Overall heat losses from the total system ranged from 250 to 1500 watts, 
depending on the operating temperature. 

3.1.3 Heat Pipe Description - Definition of Burnout 

The heat pipes tested in this program (Figure 3.1.3 and 3.1.4) used three dif- 
ferent evaporative surfaces: a single layer of 200-mesh stainless screen, 

a proprietary boiling surface — High Flux which is manufactured by Linde Division 
of Union Carbide, and a bare tube surface. The condenser surface in all cases 
was bare, as-received Type K copper tube. 

Most of the open tube thermosyphon tests described in Section 3.2.1 were per- 
formed with a bare evaporator wall of Type K copper. Although some of the data 
for benzene and P-11 were taken with the previously mentioned 200-mesh screen sur- 
face, these data have not been separately identified because no significant effect 
of the screen on critical axial heat flux was identified. Subsequent tests per- 
formed on the cocurrent flow heat pipe used either a bare evaporator tube wall or 
the Linde High Flux surface. The particular evaporative surface used in each test 
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3-8 


sequence is identified in Section 3.2.2. The Linde High Flux surface was a 
sintered layer of copper on a substrate of cupronickel. 

In all tests, the liquid charge used was one which would fill approximately 
60% of the evaporator under static saturated conditions at 25°C. Where 
relevant, actual gram loadings are called out in subsequent sections. 

Heat pipe "burnout" was defined as follows. It was observed in all cases that 
a plot of the overall heat pipe conductance versus heat flux showed a steadily 
increasing value of conductance with heat throughput until burnout was imminent. 
Near burnout, the heat pipe conductance would stabilize and then begin to 
fall, in many cases, very sharply. For all test runs, burnout was identified 
as the heat flux at which the overall conductance peaked. This means that at 
"burnout" the heat pipe was still functional; in many runs the difference 
between this burnout flux and a thermal runaway heat flux was quite small — in 
other cases, the flux at thermal runaway was as much as 50% greater. 

The heat flux at peak conductance was defined as burnout because it was felt 
that the subsequent reduction in conductance with heat input was the most clear 
signal of a fundamental and degrading change in heat pipe operation. 


3.2 Experimental Performance 

An extensive series of tests has been performed to provide an experimental basis 
for the heat pipe model developed in this project. These tests also provided a 
comparison between open core heat pipes and the enhanced capacity cocurrent heat 
pipes which were developed as part of this project. 

3.2.1 Open Core Heat Pipes 

Figure 3.2.1 is a plot of critical heat flux versus tilt from horizontal, for 
four different fluids, in a 1.91 cm I.D. bare tube heat pipe. Table 3.2.1 
presents physical properties cf the test fluids. 



3-9 



ANGLE FROM HORIZONTAL 


R-I2 
R.I14 
R-l 1 
benzene 


( 3 . 2 . 1 ) 



ANGLE FROM HORIZONTAL 


CRITICAL AXIAL HEAT FLUXES FOR A 1.91 CM I.D. 
REFLUX HEAT PIPE WITH FOUR DIFFERENT WORKING 
FLUIDS (FIGURE 3.2.1) AND CRITICAL AXIAL HEAT 
FLUX FOR A 2.54 CM I.D. REFLUX HEAT PIPE CHARG- 
ED WITH R-ll (FIGURE 3.2.2). 
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TABLE 3.2.1. TEST FLUID PHYSICAL PROPERTIES 


Working Fluid/ 
Temperature 

Liquid 
Density, p 
(g/cm 3 ) 

Vapor 
Density, P 
(g/cm 3 ) 

Bond 

Number^ 2 ) 

Vapor Energy 
Density, p h, 
W/c.3,’- c * 

Bcnzene/25°C 

0.868 

0.000505 

10.5 

0.218 

R-11/30°C 

1.46 

0.00717 

17.2 

1.26 

R-114/35°C 

1.43 

0.0213 

21.4 

2.64 

R-12/ 40 °C 

1.255 

0.0552 

24.5 

7.09 


(1) Temperature at peak heat flux used for property evaluations. 

(2) For a diameter of 1.905 cm. 


The tests clearly show four qualitative features of reflux heat pipe operation 
and the effect of fluid selection: (1) a high latent heat and gas phase density 

result in greater transport capacity at all tilts; (2) heat flux versus tilt 
behavior is very similar for all fluids tested, indicating a common physical 
mechanism of burnout; (3) heat fluxes at near-horizontal are very low; and 
(4) the absolute maximum axial flux occurs at a reflux angle less than 90°. 


Aside from the obvious beneficial effects of reducing vapor/liquid shear through 

an increased vapor phase energy density (p h, ) , it is interesting to note the 

g f g 

strong positive influence of liquid/vapor stratification which occurs when the 
heat pipe is not in a vertical orientation. The condensate stream or puddle 
formed by tilting shields more of the condensate from vapor shear by improving 
the ratio of shear area to condensate flow cross-sectional area. 


The Bond numbers of the fluids tested range from about 10 to 25. It might be 
expected that fluids with a large Bond number would have higher peak heat 
fluxes because of the weakness of surface tension forces relative to gravity 
forces, i.e., high Bond number fluids would have puddles with a minimum vapor 
shear area. The peak heat fluxes of each fluid do, in fact, fall in the same 
oi er as their Bond numbers, and the reflux angle at the absolute maximum heat 
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flux decreases with increasing Bond number. Benzene, with a Bond number of 
about 10, has a maximum axial flux of about 750 watts at a 30 to 50° tilt, 
while V.— 12, with a Bond number of 24.5, has a maximum axial flux of about 
180C catts at a tilt angle of 10 to 20°. In all cases, the axial flux at a 
90° tilt is substantially less than the flux at the optimum angle, that is, 
as much au one- third to one-half lower than the flux at the optimum tilt. 

There arc no analytical or empirical models in the literature that predict 
♦■his behavior. 

The peak heat flux in a vertical orientation has been studied in some detail, 
as discussed in Section 2.1. Figure 2.1.3 compares the data shown in Figure 
3.2.1 with the model of Tien and Chung, and the recommended correlation, equa- 
tion (2.1-10), given in Section 2. In general, the data fall within 50% of the 
Tien and Chung model, and is in good agreement with equation (2.1-10) if C^ is 
in the range 0.725 sC 0.95. 


Figure 3.2.2 gives the experimental results for critical heat flux for R-ll in 
a 2.54 cm I.D. heat pipe incorporating the High Flux evaporative surface. Also 
shown in this figure is the value for a 90° tilt predicted by equation (2.1-10). 

3.2.2 Enhanced Performance Heat Pipes 
Capacity Enhancement 

It is apparent that one of the primary factors limiting peak heat transfer in 
a reflux heat pipe is countercurrent vapor/liquid shear. The shear process 
and related phenomena of wave action and droplet entrainment serve to starve 
the evaporator of liquid and cause burnout. Hence, the primary goal must be 
to develop a means of counteracting the interfacial shear or actually using 
it to enhance performance. 
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The roost obvious means of enhancing capacity is to separate the liquid and 

vapor flows so that the destructive mechanisms related to shear are eliminated. 

For nominally horizontal reflux heat pipes, this can be done by placing a 

smaller diameter tube in the heat pipe and filling the heat pipe to a level 

such that the condenser end of the heat pipe is flooded, thereby providing a 

( 12 ) 

fluid return to the extreme evaporator end. A second technique that has 

been used on vertical reflux liquid metal heat pipes is to provide axial flow 

channels or "gutters" for condensate to minimize the vapor/liquid contact area 

(32) 

and thereby limit shear. Both of these techniques will provide substantial 

capacity improvements, and the secondary tube approach is reasonably cost- 
effective and immune to boiling within the condensate return. The axial flow 
channel is satisfactory for liquid-metal heat pipes but could pose problems 
for organic heat recovery heat pipes where the wall superheat is greater, the 
tendency towards nucleate boiling is mote pronounced, and liquid surface ten- 
sions lower. Section 3.2 describes the class of flew separator approaches 
tested in this program. These tests resulted in a separator optimized for 
organic and fluorocarbon working fluid heat pipes. 

A summary of test results for various vapor-liquid separators, along with data 
for pipes without separators, is given in Table 3.2.2. 

General Description of Test Phases 

All tests listed in the table were performed with R-ll working fluid, and are 
identified by axial transport capacity and the corresponding vepor core tempera- 
ture. Test sequences 1 through 5 were performed with a 1.905 cm I.D. bare 
evaporator surface and a relatively small 0.45 cm I.D. condensate return tube. 

The Phase 1-3 tests used a relatively short separator plate that penetrated 
only a short distance into the condenser zone. These separator plates increased 
capacity at low tilt angles by about a factor of two relative to an open- 
sectioned heat pipe, but had no substantive effect on vertical reflux performance. 



TABLE 3.2.2. DATA SUMMARY FOR CONVENTIONAL AND HIGH- PERFORM* 
REFLUX HEAT PIPE TESTS USING R-ll 1 


Tilt 

Angle 

(deq) 

Screened 

Evaporator 

(W/°C) 

Base 

Tube 

(W/°C) 

Phase 1 
Tests 
(W/°C) 

Phase 2 
Tests 
(W/°C) 

0 


<200/19 



5 


610/21 3 

>800/34 


10 

905/27 

990/28 

1100/31 

1560/43 

20 

>940/27 

1030/27 

1290/36 

>1980/50 

30 

1075/30 

1000/30 

1870/41 


40 

>1035/30 




45 


1000/30 

1710/41 

2170/48 

50 

1020/31.7 




60 

1050/30.5 

840/28 

1760/42 


70 

790/27 

831/24 

1760/41 


80 

800/29 

660/22 

1450/40 


90 

810/26 

750/23 

1340/38 

1230/30 


Phase 3 
Tests 
(W/°C) 

Phase 4 
Tests 
(W/°C) 

Phase 5 
Tests 
(W/°C) 

Phase 6 2 
Tests 
(W/°C) 


1210/36 

970/37 

No 



1880/47.5 

1730/48 

change 

3000/40 

'v-6! 

2510/56 

2400/54 

from 

3500/38 

I 

j 

>2790/61 

2840/58 

Phase 



>2720 


4 


j 


>2980/57 



>i<M 

>2720 





2710/58 

>2980/57 



] 

2440/52 

>2980/56 




1250/34 

>3030/57 


1550/32 

>io, i 


‘initial liquid load - 240 g of R-ll. 
2 Initial liquid load = 450 g of R-ll. 
3 4 . 3° tilt. 




Phase 9 2 
Tests 
(W/°C) 

^5600/75 5500/75 

>7300/82 7100/81 

>8500/81 
>14,000/74 
>9800/80 
>13,800/70 

>15,200/70 

>14,200/67 



FOLDOUT FRAME 
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TABLE 3.2.2. CONTINUED 

Identification of Tests 

Base Tube: Bare, 1.9 cm inside diameter heat pipe. Figure 3.1.3. 

Screened Evaporator: 1.9 cm diameter with 200-mesh screen lining in evaporator. 

Phase 1: 1.9 cm heat pipe with 17.9 cm long separator tube and 

0.45 cm I.D. return tube. Figure 3.1.3. 

Phase 2: 1.9 cm heat pipe with 22.8 cm long separator plate and 

0.45 cm I.D. return tube. Figure 3.1.3. 

Phase 3: 1.9 cm heat pipe with transparent adiabatic section. 

Phase 4: 1.9 cm heat pipe with transparent adiabatic section; 

105 cm long separator plate and 0.45 cm I.D. return tube. 
Figure 3.2.4. (Phase 4 modification is shown with Phase 5 
evaporator spiral. Phase 4 condensate tube is the same as 
that used in Phase 4, minus the spiral.) Figure 3.1.3. 

Phase 5: 1.9 cm heat pipe with transparent adiabatic section, 

105 cm long separator plate and spiral insert in the 
evaporator section. Phase 4 modification. Figure 3.1.4. 

Phase 6: 2.54 cm heat pipe with Linde High Flux* evaporator surface 

and no separator. Figure 3.1.4. 

2.54 cm heat pipe with High Flux evaporator, 108.6 cm long 
separator plate and 0.9 cm I.D. return tube. Figure 3.1.4. 

2.54 cm heat pipe with Hgh Flux evaporator, 108.6 cm long 
separator, oriented vertically and 0.9 cm I.D. return tube. 

Phase 9: 2.54 cm heat pipe with High Flux evaporator, 108.6 cm long 

separator and perforated 0.9 cm I.D. return tube. 

Phase 10: 2.54 cm heat pipe with High Flux evaporator, 108.6 cm long 

separator and 1.2 cm I.D. return tube. Figure 3.1.4. 


* Commercial designation of proprietary porous boiling surface manufactured by 
Linde Division of Union Carbide Corporation. 


Phase 7: 
Phase 8: 
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FIGURE 3.2.3 

Separator test sequence No. 3, 50° tilt, 2400 watts, 1/500-second exposure 
time. Condensate capture pocket on right side of separator plate is full 
of liquid, and vapor-liquid flow in heat pipe core is a mixture of annular 
and puddle flow. Heat pipe core flew patterns at 30 and 40° tilts look 
similar to the flow pattern in this photograph. 
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The Phase 3-5 test sequences were performed with a glass adiabatic section so 
as to observe heat pipe dynamics. Figure 3.2.3 is a representative photograph 
of the Phase 7 heat pipe at a flux of about 2500 watts. 

In Phase 4, the separator plate length was increased so that almost the entire 
condenser was split into halves. This permutation had a dramatic effect on 
vertical operation. With the long separator plate, vertical reflux heat pipe 
capacity exceeded the capabilities of the pumped water loop, providing a factor 
of 2.5 gain over the earlier open tube tests. The longer separator 
plate reduced transport capacity at tilt angles less than or equal to 30° by 
about 5 to 20%. Overall, the gains associated with the long separator plate 
far outweighed the small reduction in capacity at low tilt angles. In Phase 5, 
a spiral insert was used in the evaporator in an attempt to improve evaporative 
heat transfer coefficients. The spiral insert configuration is detailed in 
Figure 3.1.4. There was no measurable improvement in evaporator conductance. 

In Phases 6 through 10, a 2.54 cm I.D. heat pipe was used with an evaporator 
that was coated with the Linde High Flux porous boiling surface. Phase 6 tests 
were performed with an open core (that is, no separator plate or condensate 
return tube) to establish a performance baseline. 

Phase 7 and Phase 8 tests were performed using a 0.9 cm I.D. condensate return 
tube and a full-length separator pj.ate oriented horizontally and vertically in 
Fhases 7 and 8, respectively. No significant difference in axial transport 
capacity was found at low tilts; the orientation of the separator plate was of 
secondary importance. 

In Phase 9 tests, the 0.9 cm I.D. condensate return tube was perforated 
periodically in the evaporator zone to allow liquid to bleed onto the evaporator 
surface in a roughly uniform fashion. This should result in a capacity increase 
since all of the condensate does not have to travel to the far end of the evapor- 
ator, resulting in a reduction in liquid pressure drop. A capacity increase 
of up to 60% was observed at low tilt angles. 
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In Phase 10, a plain tube of 1.2 cm I.D. was used as the condensate return. 

The transport capacity of this heat pipe was very close to the capability 
of the Phase 9 heat pipe. The Phase 10 heat pipe is close to having an opti- 
mum ratio of condensate tube diameter to heat pipe diameter. In vertical 
reflux operation, the Phase 10 heat pipe demonstrated a transport capacity 
that is a factor of 9 greater than possible with an open cross-section heat 
pipe, while at a tilt angle of 10°, a performance increase of 2.3 times was 
obtained . 

In the following text and figures, the more significant test sequences will 
be discussed. 

Detailed Discussions 

Experimental values for burnout power, or critical heat flux, versus tilt 
angle are plotted in Figures 3.2.4 and 3.2.5 for R-ll and benzene in the Phase 10 
heat pipes. Also shown on these plots are the predicted results based on the 
Section 2.2 modeling. The results are in excellent agreement considering the 
difficulty of measuring the exact value of burnout power. In Figure 3.2.5, 
two of the data points are indicated with "up" arrows . The data point at an 
80° tilt was limited by the open pumped water loop temperature approaching 
100°C, while the point at a 90° vertical orientation is the lower limit heat 
flux obtained during the oscillations discussed subsequently. The three "up" 
arrows shown in Figure 3.2.4 were caused by limitations in the power available 
for heating the evaporator. 

In Figure 3.2.6, all burnout data for benzene and R-ll that was taken with the 
Phase 3, 4, 7, and 10 separators is compared as a function of tilt angle 
with the theoretical predictions. The experimental data a: e predicted to 
within ±25%. This level of agreement is quite satisfactory considering the 
accuracy limitation of the Lockhart-Martinelli and Owens models of two-phase 
pressure drop, and the uncertainty in identifying burnout heat flux. Exit and 
entrance qualities for the evaporator and condenser were set at 0.9 and 0.67, 
respectively, for all model calculations. 
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CRITICAL AXIAL HEAT FLUXES FOR PHASE 10 COCURRENT HEAT 
PIPE CHARGED WITH R-ll (FIGURE 3.2.4), AND WITH BENZENE 
(FIGURE 3.2.5). 






ANGLE FROM HORIZONTAL 


FIGURE 3.2.6. A COMPARISON OF EXPERIMENTAL CRITICAL HEAT FLUXES 
WITH CALCULATED VALUES FOR THE PHASE 3, 4, 7, AND 
10 COCURRENT HEAT PIPES USING BENZENE AND R-ll 
WORKING FLUIDS. 
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When operated in a vertical or near-vertical orientation, the highest capacity 
evaporator system. Phase 10, tended to oscillate between two different modes of 
operation. This oscillation is shown in Figure 3.2.7, where the evaporator 
jacket temperature is plotted versus time. The period of the oscillation is 
approximately 24 minutes in this case. A check was made of the time constant 
associated with heating and cooling the entire heat pipe heating loop system 
and it was determined both experimentally and theoretically to be approximately 
4 minutes. Thus, it appears that the observed oscillations were fundamental 
to the operation of the heat pipe and were not an artifact of the test setup. 

The experimental data taken with Phase 10 indicate that the oscillations in heat 
pipe performance were due to changes in the thermal conductance of the condenser 
section. Bezrodny (17) has observed similar oscillations in the operation of 
bare tube thermosyphons . He attributes the oscillations to changes in the 
nature of the flow in the evaporator section. It was not possible in these 
tests to establish the cause of the oscillations in the Phase 10 he-' pipe, and 
it is, therefore, not clear whether the qualitative observations of Bezrodny 
explain the changes noted here. 

Figure 3.2.8 is a plot of overall conductance versus power for the 2.54 cm heat 
pipe using the Phase 10 separator. The working fluid is benzene. This graph, 
which is generally typical for all fluids, has several features of interest. 

First, it can be seen that the conductance versus power curves follow a nearly 
common curve until the burnout power for the particular angle is approached. 

The graph also shows that the burnout power is not sharply defined. Another 
interesting result is that the burnout power appears as a reduction in the 
conductance of the evaporator section of the heat pipe, a phenomenon also observed 
in the open core heat pipes of Section 3.2.1. 

Figure 3.2.9 shows the evaporator and condenser conductances as a function of 
heat flux. Up to the point of burnout, both the evaporator and condenser 
conductances increase with heat flux. The peak evaporator conductances of 
500 to 700 watts/°C obtained with the High Flux surface correspond to a heat 
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FIGURE i.2.7. OSCILLATORY BEHAVIOR OF HEAT PIPE EVAPORATOR JACKET TEMPERATURE 
FOR THE PHASE 10 HEAT PIPE IN VERTICAL OPERATION WITH R-ll 
WORKING FLUID. 
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FIGURE 3.2.8. PHASE 10 COCURRENT HEAT PIPE CONDUCTANCE AS A FUNCTION OF TILT AND 
HEAT FLUX WITH BENZENE WORKING FLUID. TEMPERATURE RANGE = 40-80°C. 
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2 2 

transfer coefficient of 0.59 to 0.83 W/cm K (1040 to 1460 Btu/hr-ft °F) over a 

2 2 

radial heat flux density range of 3 to 18 W/cm (9500 to 57,000 Btu/hr-ft ). 

2 

The condenser heat transfer coefficients fall in the range of 0.09-0.25 W/cm K 

2 

for a radial heat flux density range of 2.70 to 16 W/cm K. It is clear t*at 
the condenser coefficient associated with the bare tube I.D. limits the overall 
heat pipe conductance, and is an area requiring further development. 
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4.0 SUMMARY AND CONCLUSIONS 

By the introduction of a plate-and-tube separator assembly into a heat pipe 
vapor core, it has been demonstrated that axial transport capacity in reflux 
mode can be improved by up to a factor of 10. This improvement is largely 
the result of eliminating the countercurrent shear that commonly limits re- 
flux heat pipe axial capacity. With benzene, axial heat fluxes up to 1800 
2 

W/cm were obtained in the temperature range 40 to 80°C, while heat flux den- 

2 

sities up to 3000 W/cm were obtained with R-ll over the temperature range 40 
to 80°C. These very high axial capacities compare favorably with liquid 

2 

metal limits; the sonic limit for liquid sodium, for example, is 3000 W/cm 
at 675 8 C. Computational models developed for these cocurrent flow heat pipes 
agreed with experimental data within ±25%, and a computer program HPIPE was 
written to facilitate optimization and calculation of heat transport limits for 
diverse cocurrent flowheat pipe designs. 

The high axial transport limits obt able under cocurrent flow conditions can 
lead to substantial reductions in heat exchanger cost through a reduction in 
the number of heat pipes required. However, to gain full benefit of the 
increase in axial capacity, the radial heat transfer resistances associated 
with vaporization and condensation must not account for am appreciable frac- 
tion of the hot stream/ cold stream temperature difference. For the develop- 
mental heat pipe tested in this program, the evaporator conductance was 
improved through use of the Linde Highflux porous boiling surface, while no 
condenser improvement was attempted beyond the smooth-walled tube used in 
these tests. Heat transfer coefficients for the evaporator and condenser fell 
in the range 0.59*0.83 W/cm^K and 0.09-0.25 W/cm^K, respectively. 

Using these heat transfer coefficients and radial heat flux densities in the 
2 

2 to 20 W/cm range (as typical of many heat recovery applications) , it cam 
be seen that overall temperature losses for the evaporator and condenser 
could ramge from 2.5-J4°C to 22-80°C, respectively. Since most heat recovery 
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applications have overall temperature differences of 250°C or less, it is 
clear that from 10% to one-third or more of the available temperature differ- 
ence may be lost in the heat pipe. It is therefore clear that improvements 
in evaporating and condensing heat transfer coefficients will be needed in 
many cases to take full advantage of the cocurrent flow heat pipe's high flux 
capabilities. 

These heat pipes are ideally suited for heat recovery applications with large 
hot stream/ cold stream temperature differences and where radial and axial 
heat fluxes are sizable, as for example in stack gas recuperators and heat 
recovery boilers, or where relatively long heat pipes are required or indi- 
cated on the basis of performance and cost studies. 
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APPENDIX A 

COMPUTER PROGRAM "HPIPE" 
PROGRAMMED IN 

EXTENDED BASIC ?0R THE PDP-11 



A-l 


£0 DEFFN H 1 T"' =£33 . 147+. 02 089 1 8+T— , 669S U8E—3+T ♦T 
3: Ci DEFFN DO -: T> = 1 . 9052 1-. 471 017E-3*T-. 324476E-5+T+T 
40 DEFFN D2 <T> =756. 4^EXP <-£986. S.- T'* ,"T++. 300665 
5 0 DEFFN V 0 CD =2 . 5 07E-4*EXP 034 0 . 1 V 2 T> 

60 DEFFN V2 CT> =-3. 0672E-3+2. 66£48E-5^T-7. 43805E-8V|>T+6. 97587E-1 

70 DEFFN F9 <R> 

71 IF F> £ 0 0 0 «3D TO 74 

72 FNF9= 1 6-P 

73 GO TO 75 

74 FNF9= . 219622--R++. 43 — . c‘61£97E—3^L0b +5. 95156E— 3 

75 FNEND 

77 DEFFN Y9 <X > 

73 IF KG < £ ij 0 0 GO TO 81 

79 FNY9 = 1 1 +:«>♦ . 1 8 03> ♦♦ 1 . 017 05* < 1 +X++ . 25> ♦♦3 . 1 57 1 £ 

8 0 bO TO 8 8 

8 1 FNY9 = '1 + 1.1 25+X++. 75> ♦♦. £ 1 3£ 1 5+ < 1 +X++ . £5 > ♦♦£ . 3 1 326 

82 FNEND 

?U DEFFN P <T> =4. 46U8E12+EXF' <— >: Jl£. 79 ' T.' 1 ■• T++. 727083 
1 0 0 DEFFN P 9 < 2 . IJ 1 > l j £> =- p 1 ♦ < 2 ♦ S 3 ♦Cl-fl •••- B > < IJ £ - U 1 > - H ♦ C ■ + A ♦ 2 ♦ S 3 > E 

110 READ 20 ! CONDENSATE RETURN TUBE LENGTH 

111 DATA 99.7 

120 READ A 0 . P 0 'RETURN TUBE C. S. AREA AND WETTED F'ERIM. 

!£1 DATA 1.119*3.75 

130 READ 21 ! EVAPORATOR LENGTH 

131 DATA 103.2 

140 READ Hi* PI ! EVAP. C.S. AREA AND WETTED PER IN. 

141 DATA 4.107.12.21 

150 READ 22 ! ADIAEATIT: LENGTH 

151 DATA 14.0 

160 READ A2.P2 !ADIA. C.S. AREA AND NETTED REPIN. 

161 DATA 4.107*12.21 

170 READ 23 .'CONDENSER LENGTH 

171 DATA 115.6 

130 READ A3.P3 SCOND. C.S. AREA AND NETTED PER IN. 

131 DATA 5. 065. 7. 98 

190 READ 24 'SEPARATOR LENGTH IN AD I A. 

191 DATA 14 

200 READ A4.P4 'SEP. C.S. APFA AND WETTED PEPIN. IN ADI A. ZONE 

201 DATA £.405.6.43 

205 READ 25 'SEPARATOR LENGTH IN COND. 

206 DATA 94.6 

210 READ A5.P5 'C.S. AREA AND RERIM. ON EACH SIDE OF SEP. IN COND. ZONE 

211 DATA £.405.6.43 

215 READ X0.X1 'EVAP. INLET AND OUTLET QUALITY 

216 DATA O'. 999 

220 READ X2.X5 'COND. INLET AND OUTLET QUALITY 

221 DATA .5«0 

£20 READ T 2 ! INITIAL SATD. VAP. TEMP. 

321 DATA 80 

240 READ H1.H2 'MIN. AND MAX. HEAT FLUXES ■ N ■ 

£41 DATA 1 ij 0 . 2 5 0 0 0 

£5 0 READ HO 'A-AIL. LIQUID *~EAB '-'HEN OPERATED VERTICALLY 
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£51 BATA 99.7 

£60 READ AS STILT ANGLE W.R.T. HOF'IZ.* DEGREES 
£61 DATA 90 

300 READ G ! GRAV. ACC. 

301 DATA 980.6 

3 1 0 A9=S I N 0 . 0 1 74533* AS) 

3£0 Z9=H0*A9 'AVAILABLE HEAD 
33 0 T=T£+£73 .15: IJ7= . 0 0 1 : H9»FNH <T> 

34 0 L 1 =H 1 ■••• H9 : L£=H£ - H9 

349 PRINT 

350 PRINT" Q 01.1) MASS FLOW 

351 PR I NT " BP-EVRP DP-ADIA. 

35 £ PRINT" 


PRESS. HEAD 
DP-SEP.-ADIA 


TILT 

DP-TIJBE 


TEMP. " 
EP-CONB' 


41 

0 

V 0 

=FNV 0 <T 

:■ : Y£=FNV£ 

T ) : l.l 0 

= 

FNDO OT 

) : l,l £ 

=FND£ 

:T) 





41 

cr 

A=l 

1 0-l>l£— 1 

: HI =l.i.i 0*b* 

Z1*A9-- 


X 1 -X 0) 








4£ 

o 

F 0= 

= . 5*Z5 - 

Z3: X3=X£+ 

FO* 0X5 

- 

X£) : X4 

_ y Cj __ 

0X3-X2) : F 

C = 

l.l £ ♦ 0 1 X 

1 - 1 ) 

♦ ♦£,"1,1 0 

43 

o 

X6- 

= 1 +A*X£ 

: X7= 1 +A*X 

4 : X8= 1 

- 

X 1 : X9= 

1 +A* 

XI 






44 

o 

co= 

=4 -- OPO* 

VO) : C 1 =4 

<P 1 ♦Y 0 

I 

: C£=4 

0P£* 

Y d.' ■ L •. 

= 4 / 

OP 

5*YO) :C 

4=4-' 

0F'4*Y£ 

45 

o 

r-,= 

=4 - <P£* 

Y0 ) : C6=4 •• 

0F'4*Y 0 

’l 

: C7=4- 

0P5* 

0 £) : CI: 

=4." 

OP 

1 ♦Yc' 1 



46 

0 

B0= 

= . 5*P 0* 

Z 0 ■: l>t 0 ♦A 0 

♦♦3) : E 

1 

=.5*P1 

s (lil 0 

♦A1 

> : E 

c. = 

. 5*P£ -' '■ 

l.l £ ♦ R 

£*♦3) 

47 

0 

E3= 

= . 5*P5-- 

0W£*A5**3 

> :B4=. 

5 

♦F'4*Z4 

0 1 i.l 3 

♦A4**0 

> : E 

CJ = 

El* *1.1 0 ••• 

l'l£) 


47 

4 

E4 = 

=X9 <£♦ 

,ilO*A4*A4) 

: E5=E4 

♦ 

0A4--A5 

:■ ♦♦£ 







47 

c 

IF 

Z£ <' . 0 1 

bD TD 49 

C 











476 Z6-Zc~Z4 

477 IF Z4 < . 01 bO TD 490 

478 IF ABS <Z4-Z£> <.01 GO TD 485 
430 E£=E4: E=R4 •A£: GO TD 500 
435 E£=E4 : E=A4 H 1 : bD TO 500 
490 E£=E5: E=A5-- A£: bD TD 5 00 
495 E£=E5: B=A5 • A1 : Z6=0 

500 E£ = . 45* l-B) *E£: E3=30* o £+A* 0X3+X4) > <W0*A5*A5> ! ENTRANCE EXIT FACTORS 
510 b 1 = - 1 0 l.l £ ♦ F N P •: T ) ♦ A 1 ♦ A 1 > : V 1 = 1 0 l.l 0 ♦ A 1 ♦ A 1 > 

515 IF Z5<. 01 bD TD 550 

53 0 F5=iii£* 1 1 X£- 1 > ♦♦£.-•• l.l 0 : F6=W£* 0 1 X4- 1 > ♦♦£,•• t.i 0 
5£5 b5=b 1 ♦ o A 1 ■ A5) ♦♦£ : V5= 1 •••" Nil 0*A5*A5> 

5 3 0 H 5 = U 0 ♦ b ♦ Z 5 ♦ A 9 s 0 X 3 — X £ ) : H 6 = H 5 ♦ 0 X 3 — X £ ) s' 0 X 5 — X 4 1 

54 0 0 5 = L . 0 b <1+A*X£> 0 1 +A*X 3 > > : 0 . 6 = L 0 b < 0 1 + A ♦ X 4 ) < 1 + A ♦ X 5 > > 

55 0 0 1 = L D G o l + A ♦ X 0 > 0 1 + A ♦ X 1 > > 

56 0 B=A0/A5s E0=. ££5* 0 1 -B> <M0*A0*A0> : E 1 = . 45*B**3-" 0 1 +B> 

57 0 E 1 = 1 E 1 + 1 < 1 +B) ♦♦£> £ ♦ I « 1 0*A 0*A 0) : E6=6 0*E*B - <l.«» 0*A 0*A0* 0 1 +B> ♦♦£> 

590 U6=- 1 : U9= 1 E- 1 £ 

6 0 0 M= 0 L 1 +L£) £ s M4= •: L£-L 1 > s4 
610 FDR 1=1 TD 100 

6 1 5 M 0=M*X8 : M 1 =M*X 1 : M£=M*M : M6=F 0*M 1 : M8=M 1 ♦Hi 

7 0 0 ! CONDENSATE TUBE PRESSURE DROP- FRICT. AND ENTR. 

710 F=FNF9 C 0*M> 

7 £ 0 S 0 = B 0 ♦ M £ ♦ F : S 1 = • E 0 + E 6 ♦ K N F 9 e 4 ♦ M • 0 V ' 0 ♦ < P 1 - p 0 > > • + E 1 > ♦ M £ 

300! ADIABATIC ZONE 
3 05 IF Z£0.Ol bD TD 395 

8 1 0 F = F N F 9 C £ ♦ M 1 ) : X = F £ ♦ F N F 9 1 C 5 ♦ M 0) ■■■" F : !■ 8 = C 5 ♦ M 0 : V = F N V 9 < X ) 
d 1 ' 6 3 * Y ♦ B c! ♦ M 8 ♦ F 
: '■ 0 “5= 63*Z6: GO TC 900 
; ^ ‘ 5 = 0 



9 0 0 ! E V fiP OP ATOP ZONE 

905 F*FNF9 CC8+M 1 1 : X*F£+FNF9 « C 1 ♦MO ■••• F : K 8=C 1 ♦M 0 : Y=FNY9 <X> 

9 1 0 B“G 1 ♦ME : C=V 1 ♦MS : B6= 1 +B^X 0 : B7=l +B*X 1 

911 IF B6<0 GO TO 1 £9 O • IF B7<0 GO TO 1 £9 U 
915 R 1 = L 0 G < E 7 •••■' E 6 > 

9£0 S 3 = - V ♦ B 5 ♦ M 8 ♦ F X 9 ! INLET BP-DZ COMPAT. WITH EXIT BP-DZ 

93 0 $4=FNP9 <Z 1 > X 0 » X 1 > 

94 0 S£=H 1 ♦ < i? 1 +R 1 > <B-A> ! GRftV . BP 

1 0 0 0 ! SEPARRTOR-AB I fl . SECT I ON 

1004 S7=E£^M£ ! ENTRANCE 

1005 IF Z4C.01 GO TO 1095 

1 0 l 0 F=FNF9 <C4#M 1 > : X=F£^FNF9 C6^M 0> • - F : K8=C6^M 0 : Y=FNY9 OO 
1030 S6*Y^B4^H8^F ! FRICT. DP 

1090 GD TD 1100 


1095 36=0 

1 1 0 0 ! CONDENSER ZONE 
1105 IF Z5 < . U 1 GO TO 1195 

1110 M 3 = M 1 X £ : M 5 = M 3 ♦ Cl - X £ ;• : M £ = M 3 ♦ M 3 

1111 P 4 = E 3 ♦ M £ ♦ F N F 9 C C 3 ♦ M 3 > ! EXIT ENTRANCE RT PLATE TERMINATION 

1115 B=G5*M£ : C=V5+M£ : B6= 1 +B^X£ : B7= 1 +B+X3 : B3= 1 +E^X4 s B9= 1 + B+X5 

1116 IF B6\0 GO TO 1£90* IF B7<0 GO TO 1£90 

1117 IF B3 < 0 GO TO 1£90: IF B9 0 GO TO 1£90 

1118 R 1 = LOG ■: B7--'B6> 

1 1 £ 0 F=FNF9 <C7+M1> : X=F5+FNF9 <C3+M5> - F : K 8 = C 3 ♦ M 5 ' V = F N Y 9 <X> 

1 1 £5 S 3 = - Y ♦ B 3 ♦ M 8 ♦ F X 6 : P5=FNP9 •: Z5 - X£ . X3J 
1 1 3 0 R6=H5+ < P 1 +Q5) 1 B— A ! GF’AV . DP 
115 0 M5=M3+ •: 1 -X4> : R 1 =LOG ■ B9 B8> 

1 1 o u F=F N F 9 1 C 7 ♦ X 4 ♦ M 3 ' • X = F 6 ♦ F N F 9 C L 3 ♦ M5 1 ’ Y = F N Y 9 1 X J 

1165 33=-Y+B3+ 1 M 3+X4) ++£+F -"X7 : R7=FHP9 , ’Z5 » X4 « X5> ! FRICT . + ACC. BP 
1 1 7 U F'8=-H6+ 1 L'6+R 1 > 1 B—A.) ! bPAV . DP 

1190 GO TO 1 £ 0 0 
1 1 95 P5= 0 ; P6= 0 : P7= 0 : R8= 0 
1 £ 0 0 ! P P E S S l J RE BALA N C E 

1310 33=1.1 0*G»Z£^A9-‘X9 ! GRAY' BP IN ADI A. 

1££0 39=1.1 0*G^Z9 


:n N = \’-i- 


34-3 £-3 1 -3 0 


1 £ 4 ij j = \ — R 4 — R 5 ~ R 6 — R 7 — R 8 
1 c 5 0 IF v < i.i h O TO 1 c' 9 0 

1 £6 0 M=M+M4 : IJ8=AB 3 <: M - U9 1 - 1 
1 £ 7 0 IF A B 3 I j 3 > ■ IJ 7 GO TO £ 0 0 0 
1 £ 8 0 I J 9 = M : G O TO 13 0 0 
i £9 0 N=M-M4: U6=l 
1300 M4=M4 . £ : NEXT I 

£000 SPRINT STATEMENTS 


£001 IF IJ9 1.1 0E- 1 £ GD TO 3 0 00 

£ O O £ IF U 6 < U h o r O J 0 1 0 

£010 C=M 1 ♦H9 'LATENT HEAT TRANSFER 

£ ij £ IJ P R I N T i. 1 * M * j 9< Ay « T £ 

£ 0£5 IJ 1=3 £+3 4 : U£= 3 5! IJ3=3 6+3 7 ! 1 J 4=3 0+3 1 ’ IJ5=P4+R5+P6+P7 +R'3 
£ 0 30 PRINT IJ 1 * U £ , IJ 3 , IJ 4 « IJ 5 
£050 PRINT 
£099 GO TO 9999 

3 0 00 PP I NT "MAS I FLOW LESS than min. VALUE"! GO TO 9999 
3 •■■10 PR I NT " MAS 3 FLOW EXCEEDS MAX. VALUE" 



